NASA 

Technical  Memorandum  4186 


AVSCOM 
Technical  Report  88-A-003 


Rotor  Systems  Research  Aircraft 
Risk-Reduction  Shake  Test 


J.  Brent  Wellman 
Aeroflightdynamics  Directorate 
USAART A- AVSCOM 
Ames  Research  Center 
Moffett  Field,  California 


NASA 

National  Aeronautics  and 
Space  Administration 

Office  of  Management 

Scientific  and  Techmca, 
Information  Division 


'■  y-rin  a 


'7 


1990 


SUMMARY 

s' 

A  shake  test  and  an  extensive  analysis  of  results  were 
performed  to  evaluate  the  possibility  of  and  the  method 
for  dynamically  calibrating  the  Rotor  Systems  Research 
Aircraft  (RSRA)!  The  RSRA  airframe  was  subjected  to 
known  vibratory  loads  in  several  degrees  of  freedom  and 
the  responses  of  many  aircraft  transducers  were  recorded. 
Analysis  of  the  transducer  responses  using  the  technique 
of  Dynamic  Force  Determination  showed  that  the  RSRA, 
when  used  as  a  dynamic  measurement  system,  could  pre¬ 
dict,  a  posteriori,  an  excitation  force  in  a  single  axis  to  an 
accuracy  of  about  5%  and  sometimes  better.  As  the  anal¬ 
ysis  was  broadened  to  include  multiple  degrees  of  free¬ 
dom  for  the  excitation  force,  the  predictive  ability  of  the 
measurement  system  degraded  to  about  20%,  with  the 
error  occasionally  reaching  100%.  The  poor  performance 
of  the  measurement  system\is  explained  by  the  nonlinear 
response  of  the  RSRA  to  vibratory  forces  and  the  inade¬ 
quacy  of  the  particular  method  used  in  accounting  for  this 
nonlinearity.  The  RSRA  shake  test  has  revealed  areas  of 
concern  for  an  RSRA  dynamic  calibration,  and  has 
pointed  out  weaknesses  in  the  force-determination 
method..  Suggestions  are  made  to  extend  the  method  to 
nonlinear  systems  both  by  using  a  force-determination 
variant  and  by  approaching  the  problem  using  a  totally 
different  model  of  the  airframe  response. 

INTRODUCTION 

The  NASA/Army  Rotor  Systems  Research  Aircraft 
(RSRAj^s  a  flying  test  bed  for  in-flight  investigation  and 


verification  of  new  helicopter  rotor  system  concepts  and 
technologies  .(fig.  1).  It  is  designed  specifically  to  provide 
research  information  similar  to  that  obtained  from  wind- 
tunnel  testing  of  rotors,  but  its  data  are  unique  because 
they  are  obtained  in  flight.  A  complete  description  of  the 
aircraft  is  given  by  Hall  and  Merrill  (ref.  1). 

The  RSRA  is  unique  among  rotorcraft  in  that  it  pos¬ 
sesses  a  rotor-loads  balance  system  consisting  of  seven 
load  cells  mounted  beneath  the  main  transmission;  the 
load  cells  are  designed  to  detect  and  measure  forces 
induced  by  the  rotor  system  on  the  rotor  shaft.  The  geom¬ 
etry  of  this  force-measuring  system  (fig.  2)  is  such  that 
vertical,  lateral,  and  longitudinal  loads  can  be  measured, 
as  well  as  moments  about  the  pitch,  roll,  and  torque  axes. 
In  addition,  the  transmission  housing,  transmission  frame, 
fuselage,  and  wings  are  equipped  with  accelerometers, 
and  the  main-rotor  shaft  is  equipped  with  strain  gauges. 

One  of  the  primary  research  goals  of  the  RSRA  is  to 
measure  static  and  dynamic  loads  imposed  by  rotor  sys¬ 
tems  on  the  RSRA  rotor  shaft.  Such  measurements  will 
provide  the  rotorcraft  community  with  a  national  facility 
for  testing  new  rotors  and  rotor  systems  concepts.  The 
static  calibration  of  the  RSRA  balance  system  is  described 
by  Acree  (ref.  2)  and  will  not  be  treated  here.  The  method 
of  a  dynamic  calibration  of  RSRA  was  the  subject  of  a 
study  conducted  at  NASA  Ames  Research  Center  by  the 
Rotorcraft  Flight  Investigations  Branch. 


Figure  1.-  Rotor  Systems  Research  Aircraft. 


Figure  2.-  Main-rotor  load-measurement  system. 


It  will  be  necessary  to  conduct  a  lengthy,  and  costly, 
dynamic  calibration  of  the  RSRA  loads-measurement 
system  to  provide  the  RSRA  with  the  capability  to  mea¬ 
sure  in-flight  vibratory  loads.  This  calibration  will  require 
the  application  of  a  series  of  known  vibratory  loads  on  the 
aircraft  and  the  measurement  of  the  aircraft  structural 
response  to  those  loads.  The  resulting  data  may  then  be 
correlated  statistically  to  yield  dynamic  calibration 
matrices. 

A  risk-reduction  shake  test  was  conducted  to  define 
the  validity  and  acceptability  of  the  proposed  vibratory- 
load  measurement  technique,  called  Dynamic  Force 
Determination.  This  document  presents  results  of  the  test, 
interprets  the  results  in  terms  of  RSRA  airframe 
dynamics,  and  discusses  directions  for  further  work  to 
evaluate  the  potential  for  dynamic  load  measurement  with 
the  RSRA.  The  appendix  provides  detailed  results  for  all 
conditions  of  the  test. 

SYMBOLS 

[A]t  matrix  transpose  of  A 

[A]+  matrix  pseudoinverse  of  A 

C  calibration  factor 


cn  nth  coefficient  of  a  polynomial 

Fi  force  amplitude  in  the  ith  degree  of  freedom 

F.  predicted  force  amplitude  in  the  ith  degree  of 
freedom 

Fm  complex  force  of  master  shaker 

Fs  complex  force  of  slave  shaker 

Hjj  transfer  function  of  the  jth  transducer  response  to 
the  ith  degree  of  freedom  of  the  inputs 

ij,k  indices 

Im(z)  imaginary  part  of  z 

n  number  of  test  runs  for  a  given  load  direction 
Re(z)  real  part  of  z 

Rj  response  of  the  jth  transducer 

e;  normalized  error  in  the  ith  degree  of  freedom 
\  error  radius  of  predicted  force 
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summation 


I 

0  phase  angle 

to  frequency 

BACKGROUND 

The  RSRA  force-measurement  system  consists  of 
seven  load  cells  mounted  at  the  base  of  the  transmission 
as  shown  in  figure  2.  Its  application  to  the  measurement 
of  static  loads  is  discussed  in  detail  by  Burks  (ref.  3)  and 
Acree  (ref.  4). 

For  the  purposes  of  force  measurement,  a  hub- 
centered,  nonrotating  coordinate  system  is  used  (fig.  3). 
Load-cell  outputs  may  be  transformed  into  static  loads  in 
this  frame  of  reference  by  using  a  calibration  matrix 
whose  coefficients  are  determined  by  statistical  regression 
techniques  on  calibration  data  obtained  from  ground  test¬ 
ing.  A  similar  technique  may  be  used  in  principle  to  mea¬ 
sure  simple  sinusoidal  excitation  at  a  single  frequency. 
The  method,  known  as  Dynamic  Force  Determination, 
uses  transfer  function  data  to  establish  a  complex 
calibration  matrix,  which  may  then  be  used  to  transform 
transducer  amplitude  and  phase  data  into  force  amplitudes 
in  any  given  coordinate  system.  A  more  detailed 
explanation  is  provided  in  a  later  section. 

A  number  of  extra  transducers  were  added  for  this 
test  to  allow  more  accurate  determination  of  the  forces. 
Twenty-eight  transducers  were  used  in  an  attempt  to  find 
the  best  combination  of  sensors  for  a  hybrid  force- 
measurement  system.  The  transducers  were  divided  into 
three  “sets”: 


Figure  3.-  RSRA  rotor-axis  system. 

These  sensors,  in  various  combinations,  were  evalu¬ 
ated  for  their  ability,  as  sets,  to  allow  accurate  determina¬ 
tion  of  the  applied  vibratory  load  in  five  coordinate  axes, 
as  well  as  in  a  multiple-axis  shaking  mode.  For  this  pur¬ 
pose  the  three  transducer  groups  were  combined  to  form 
five  “cases”: 

1 .  Case  1 :  Static  load  cells  alone 

2  Case  2:  Static  load  cells  plus  the  transmission 
accelerometers 

3.  Case  3:  Static  load  cells  plus  the  airframe 
accelerometers 

4.  Case  4:  Airframe  accelerometers  alone 


1.  Set  1:  Load  cells  of  the  static  force  measurement 
system 

2  Set  2:  Accelerometers  mounted  in  orthogonal 
groups  of  three  on  the  transmission  case  and  force¬ 
measuring  system  frame 

3.  Set  3:  Accelerometers  mounted  on  the  airframe 

Details  of  the  location,  type,  and  orientation  of  these 
transducers  are  given  in  table  1,  and  illustrated  in  fig¬ 
ures  2, 4,  and  5. 


5.  Case  5:  All  load  cells  and  accelerometers 

The  mathematical  techniques  presented  below  were 
used  to  evaluate  each  case  as  a  potential  dynamic  load 
measurement  system. 

Vibratory  steady-state  loads  on  a  generalized 
dynamic  system  at  any  particular  frequency  may  be 
determined  by  applying  a  linear  transformation  to  the 
complex  response,  assuming  system  linearity.  In  generic 
terms, 


Force 

input 


Some  quality 
of  the  structure 


x  Struct,,ral  1  (1) 

response  1 
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TABLE  1.- TRANSDUCER  LOCATIONS 


Mnemonic 

Type“ 

Scib 

Location 

Orientation 

MRDRAG 

L/C 

1 

Drag  load  cell 

Aft 

MRLIFTA 

L/C 

1 

Left  fwd.  load  cell 

Up 

MRLIFTB 

L/C 

1 

Right  fwd.  load  cell 

Up 

MRLIFTC 

L/C 

1 

Left  aft  load  cell 

Up 

MRLIFTD 

L/C 

1 

Right  aft  load  cell 

Up 

MRGBQCE 

L/C 

1 

Fwd.  torque  load  cell 

Left 

MRGBQCF 

L/C 

1 

Aft  torque  load  cell 

Right 

LOMGB 

Acc 

2 

Atop  main  gearbox  casing 

Fwd 

LMGB 

Acc 

2 

Atop  main  gearbox  casing 

Left 

VMGB 

Acc 

2 

Atop  main  gearbox  casing 

Up 

XMRFBPV 

Acc 

2 

Fwd  transmission  frame 

Up 

XMRFBPL 

Acc 

2 

Fwd  transmission  frame 

Left 

XMRFBPLO 

Acc 

2 

Aft  transmission  frame 

Left 

STA56NV 

Acc 

3 

Station  56,  centerline 

Up 

STA56NL 

Acc 

3 

Station  56,  centerline 

Left 

VWGTPRT 

Acc 

3 

Right  wing  tip 

Up 

VWGTPLT 

Acc 

3 

Left  wing  tip 

Up 

LIGB 

Acc 

3 

Base  of  vertical  tail 

Left 

VIGB 

Acc 

3 

Base  of  vertical  tail 

Up 

LATTPYLN 

Acc 

3 

Top  of  vertical  tail 

Left 

VTTPYLON 

Acc 

3 

Top  of  vertical  tail 

Up 

XMRFLV 

Acc 

3 

Under  left  fwd.  trans.  frm. 

Up 

XMRFRV 

Acc 

3 

Under  right  fwd.  trans.  frm 

Up 

S222FLRV 

Acc 

3 

Station  222,  floor,  right 

Up 

S222FLLV 

Acc 

3 

Station  222,  floor,  left 

Up 

S450FLV 

Acc 

3 

Station  450,  floor,  c.l. 

Up 

S450FLL 

Acc 

3 

Station  450,  floor,  c.  1. 

Left 

S450OVRL 

Acc 

3 

Station  450,  top,  c.  1. 

Left 

“L/C  =  load  cell;  Acc  =  accelerometer. 

*Set  =  number  as  described  in  text:  1  =  baseline  load  cells; 

2  =  transmission  accelerometers;  3  =  fuselage  accelerometers. 
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For  example , 


'  Re(F,) 

M  =  M  (Rj] 

(2) 

lm(F,) 

where  { F; }  is  a  vector  of  time-independent  amplitudes  of 
the  forcing  function  describing  the  structural  force  input: 

f(t)  =  Fiei(0t=Ficos(o)t  +  <j>)  (3) 

and  [C]  is  a  matrix  representation  of  the  modal  response 
model  of  the  structure.  In  this  respect,  the  structure  is 
treated  as  a  black  box  problem,  in  which  the  behavior  is 
inferred  from  input-output  measurements. 

An  approach  to  the  formation  of  the  [C]  matrix  is  to 
use  transfer  functions  of  the  transducers  as  matrix  ele¬ 
ments.  For  such  purposes,  a  transfer  function  may  be 
defined  as 

Hji =  Rj  /  Fj  (4) 

if  the  system  is  assumed  linear. 

In  the  frequency  domain,  all  quantities  in  equation  (3) 
are  complex,  so  that 

[Hji]  {Fi|  =  (Rj|  (5) 

which  indicates  the  implied  vector  nature  of  the  complex 
quantities,  and  in  which 

'RejHji)  -Im(Hjj) 

(6) 

Im(Hji)  Re(Hji) 

For  such  a  simple  one-degree-of-freedom  (DOF)  sys¬ 
tem  it  is  easy  to  see  that 

[c]=[Hji]_1  m 

If  this  idea  is  extended  to  six  degrees  of  freedom,  a  mini¬ 
mum  of  six  responses  must  be  recorded  for  each  input  to 
determine  all  six  inputs.  Equations  (S)  and  (7)  hold,  and 


Im(F6) 


a  column  vector  of  order  2  x  DOF,  and 


The  [H]  matrix  is  in  this  case  a  6  x  6  square  array  of  2  x  2 
arrays. 

Unfortunately,  equations  (2),  (5),  and  (7)  are  not  use¬ 
ful  in  an  engineering  sense  because  they  exhibit  sensitiv¬ 
ity  to  errors  in  the  measurements  of  the  elements  of  [Hji]. 
If  the  [H]  matrix  is  inverted,  along  with  any  errors  in  the 
measurement  of  its  elements,  it  may  yield  physical  infor¬ 
mation  that  is  meaningless.  The  effects  of  such  an  inver¬ 
sion  are  insidious  because  they  are  independent  of  how 
well-conditioned  [H]  is.  This  is  a  serious  problem,  but  it 
may  be  alleviated  by  the  utilization  of  redundancy,  as 
related  by  Bartlett  and  Flannelly  (ref.  5).  Redundancy  is 
achieved  by  allowing  the  row  space  of  the  response  to  be 
larger  than  the  column  space  of  the  input  forcing  function. 
In  practice  this  may  be  achieved  by  using  a  large  number 
of  transducers  relative  to  the  number  of  degrees  of  free¬ 
dom  of  the  input 

If  a  redundant  response  space  is  used,  then  [H]  is 
rectangular,  being  of  dimension  12  x  (2n),  where  n  is  the 
number  of  transducers,  and  is  in  general  singular  and  can¬ 
not  be  inverted.  To  transform  [H]  into  [C]  we  use  the 
following: 

[C]=([HJT  [H]  (10) 


which  is  the  pseudoinverse  first  introduced  by  Moore 
(ref.  6).  The  pseudoinverse  may  be  shown  to  be  a  least- 
squares  solution  to  the  redundant  transformation.  Now 
equation  (2)  becomes 
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M'KrW  <■» 

In  principle,  this  approach  should  give  the  elements  of  Fj 
to  an  arbitrarily  high  degree  of  accuracy,  assuming  linear 
responses  of  both  the  RSRA  structure  and  the  transducers. 
In  practice,  however,  the  accuracy  is  limited  by  the  large 
number  of  transducers  required  for  a  relatively  small  gain 
in  accuracy,  and  by  the  mathematical  stiffness  and  ill- 
conditioning  of  the  system  as  the  ratio  of  transducers  to 
input  DOFs  increases. 

The  concept  of  Dynamic  Force  Determination  may  be 
extended  to  cover  inputs  other  than  sinusoids  by  allowing 
the  input  force  Fj  to  be  a  function  of  co: 

{Fi(co)}  =  [Hji(co)]+{Rj(co)}  (12) 

Figure  6  shows  Hjj((o)  for  a  hypothetical  ideally  lin¬ 
ear  structure  plotted  versus  to  and  Fj.  In  practice,  if  the 
calibration  matrix  [Hjj]+  is  known  discretely  for  a  large 
number  of  frequencies  distributed  over  some  bandwidth, 
then  the  input  vector  (Fj)  may  be  input  discretely  into 
equation  (12)  over  the  frequency  band  to  reconstruct 
(Fi(co)} ,  a  spectrum  of  the  input  excitation.  The  proposed 
process  of  dynamic  calibration  is  the  determination  of  a 


suitably  large  number  of  such  discrete  calibration 
matrices,  as  determined  from  {Rj(co)}  data  taken  in  a 
shake  test,  to  cover  the  frequency  band  of  interest. 

The  actual  response  of  the  RSRA  airframe  was  not 
expected  to  be  truly  linear  because  of  the  effects  of 
mechanical  clearances  (such  as  those  around  rivets)  and  of 
breakout  friction  in  the  joints  connecting  the  transmission 
to  the  load  cells.  These  affect  the  transmission  of  vibra¬ 
tion  from  the  input  excitation  point  to  the  response  trans¬ 
ducer,  resulting  in  a  transfer  function  that  is  itself  a  func¬ 
tion  of  the  input  force.  The  effect,  called  ankylosis,  is 
such  that  Hji  is  not  constant  with  respect  to  input  force, 
but 

Hjj  =  3Rj  /  dFj  =  /(Fj)  *  Rj  /  Fj  (13) 

The  effect  was  expected  to  be  more  noticeable  at  low 
force  levels  than  at  high  force  levels.  The  expected  effect 
of  ankylosis  on  data  for  an  otherwise  linear  airframe  is 
shown  in  figure  7.  At  low  force  levels,  the  frictional 
breakout  forces  and  rivet  clearances  have  a  significant 
effect  on  the  response.  As  force  levels  are  increased, 
however,  these  effects  drop  off  relative  to  the  applied 
force  so  that  the  transducer  Rj  shows  a  net  nonlinear 
transfer  function  with  respect  to  Fj, 


Figure  6.-  Transfer  function  of  hypothetical  linear  structure,  (a)  Real;  (b)  imaginary. 
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Figure  7.-  Effect  of  ankylosis  on  transfer  function  data,  (a)  Real;  (b)  imaginary. 


At  high  force  levels,  this  source  of  nonlinearity  is 
ameliorated  and  the  response  becomes  more  nearly  linear 
with  respect  to  the  applied  force.  The  transfer  function  Hji 
therefore  approaches  some  constant  value  asymptotically 
as  force  level  is  increased.  This  asymptote  may  be  used  as 
an  estimate  of  Hji,  in  which  case  this  source  of  nonlinear¬ 
ity  is,  in  effect,  ignored.  If  ankylosis  is  indeed  small,  the 
use  of  the  asymptote  is  justified. 

To  find  asymptotes  for  each  transducer  for  each 
degree  of  freedom,  it  was  required  to  take  data  at  numer¬ 
ous  force  levels  and  to  then  find  the  best  fit  of  those  data 
to  some  asymptotic  curve.  The  curve  chosen  was 

Hji  =  c0+^  (14) 

M 

This  curve  was  chosen  for  no  other  reason  than  that  it  was 
a  curve  that  had  the  required  property  of  approaching  a 
horizontal  line  y  =  co  for  large  Fj.  The  constants  co  and  ci 
are  determined  by  least-squares  statistical  analysis  of 
bivariate  data  (i.e.,  ordered  pairs  of  force  and  transfer 
function).  Additionally,  the  horizontal  line 

Hji  =  c2  (15) 

was  used  as  an  alternate  fit.  The  constant  C2  represents  the 
best  approximation  of  the  data  that  neg’ects  the  depen¬ 
dence  of  ankylosis  on  force  level.  This  approach  consists 


of  the  taking  of  a  simple  statistical  mean  on  the  univariate 
transfer  function  data.  The  least-squares  approaches  of 
equations  (14)  and  (15)  are  shown  in  figure  8.  Equa¬ 
tion  (1 1)  may  then  be  used  with  statistically  determined  co 
or  C2  estimates  for  Hji  to  build  matrix  [C]. 


Figure  8  -  Curve-fit  approaches. 

For  the  purpose  of  this  test,  nonlinear  mass,  stiffness, 
and  damping  effects  were  not  modeled.  If  the  mass  or 
stiffness  matrices  for  the  structure  contain  elements  that 
are  functions  of  the  applied  excitation,  then  modal  fre¬ 
quencies  are  also,  in  general,  functions  of  applied  excita¬ 
tion.  The  expected  effect  of  this  phenomenon,  if  present, 
is  shown  in  figure  9.  Note  that  the  constant-frequency 
lines  in  the  vicinity  of  this  hypothetical  nonlinear  mode 
vary  with  force  in  a  manner  that  is  not  easily  modeled  by 
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Figure  9.-  Effect  of  frequency  shift  on  transfer  function,  (a)  Real;  (b)  imaginary. 


the  scheme  given  above,  whereas  constant-frequency  lines 
not  close  to  the  drifting  peak  are  relatively  well  behaved. 
Consequently,  at  frequencies  close  to  such  modes,  co  and 
C2  values,  if  indiscriminately  used  as  values  for  Hjj  in 
matrix  [C],  should  result  in  large  errors  when  applied  to 
the  input  data. 

DESCRIPTION  OF  TEST 

The  RSRA  was  tested  in  the  compound  configuration 
with  the  rudder,  lower  and  upper  horizontal  stabilizers, 
wing  and  auxiliary  engines,  and  pylons  in  place.  It  was 
prepared  for  the  test  by  removing  the  rotor  and  rotor  hub 
and  replacing  them  with  a  specially  designed  aircraft  sup¬ 
port  adapter.  The  adapter  was  designed  to  match  the  mass 
and  inertial  properties  of  the  rotor  head.  The  fuel  system 
was  drained  and  preserved,  and  the  T-58  engines  were 
replaced  with  scrap  engines  of  the  same  model.  Substi¬ 
tutes  for  the  TF-34  auxiliary  engines  were  not  available, 
necessitating  a  daily  period  of  motoring  for  each  engine, 
to  preserve  the  bearings  and  seals.  Extra  accelerometers 
were  added  to  the  transmission  and  airframe. 


The  aircraft  was  suspended  indoors  from  a  large 
movable  crane  (fig.  10).  The  suspension  system  included 
a  soft  bungee  bundle  to  provide  a  free-free  environment 
simulating  flight.  A  safety  cable  with  a  factor  of  safety  of 
6  was  used  to  back  up  the  bungees.  Hydraulic  power  was 
supplied  to  the  aircraft,  thus  allowing  the  wing  to  be 
moved  to  a  flight  attitude  and  the  landing  gear  to  be 
retracted.  The  cabin  door  had  to  be  removed  because  of 
heat  buildup  from  instrumentation,  thereby  preventing  a 
totally  faithful  re-creation  of  the  flight  configuration  of 
the  aircraft 

Excitation  was  provided  by  a  two-axis  hydraulic 
inertial  shaker  mounted  on  the  aircraft  support  adapter 
(fig.  11).  The  mount  provided  for  longitudinal,  lateral, 
vertical,  roll,  and  pitch  excitation,  as  well  as  a  combina¬ 
tion  of  lateral,  longitudinal,  and  vertical  shears  (fig.  12). 
Because  of  structural  limitations  of  the  main  transmission, 
vertical  moments  (i.e.,  torques)  were  not  applied. 
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Figure  10.-  RSRA  in  shake  test  facility. 


Figure  1 1 .-  Hydraulic  inertial  shaker  for  hub  excitation. 
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Figure  12  -  Shaker  configurations. 


The  locations  of  the  28  aircraft  transducers  used  are 
shown  in  figures  2, 4,  and  5.  The  transducers  are  given  by 
their  mnemonics  in  table  1 ,  which  shows  the  transducer 
types,  locations,  and  orientations.  The  specific  selection  of 
the  28  transducers  was  intended  to  evaluate  the  various 
proposed  approaches  for  the  measurement  of  dynamic 
forces  at  the  rotor  hub.  The  load-measurement  system 
load  cells  alone  (fig.  2)  were  chosen  to  evaluate  the  possi¬ 
bility  of  using  the  static  system  to  measure  hub  dynamic 
loads.  Two  orthogonal  groups  of  three  accelerometers 
each  (fig.  4)  were  also  added  to  evaluate  the  addition  of 
transmission  and  load-system  inertial  properties  into  the 
measurement  approach.  Airframe  accelerometers  (fig.  5) 
were  selected  to  evaluate  force  determination  without  the 
use  of  the  static  load  system  and  to  gather  modal  data  on 
the  RSRA  airframe.  Before  each  test,  the  shaker  and 
lower  support  assembly  were  moved  into  the  proper  con¬ 
figuration  for  the  particular  axis  to  be  excited.  Figure  12 
shows  the  various  configurations. 

During  each  test,  the  aircraft  was  lifted  from  the  floor 
of  the  test  bay,  and  the  landing  gear  was  retracted;  the 
wing  was  placed  at  an  angle  of  incidence  of  10°,  which  is 
typical  of  that  used  in  flight.  The  tail  landing  gear  was 
loosely  constrained  with  nylon  line  to  prevent  large 
amplitude  movements  of  the  tailcone  and  empennage 
(which  could  damage  either  the  equipment  or  the  aircraft). 

Structural  limitation  studies  were  first  made  in  each 
of  the  axes  to  determine  acceptable  test  limits  to  follow. 


These  studies  were  conducted  by  exciting  the  rotor  shaft 
at  a  low  force  level  (typically  50  lb  (222.4  N))  at  frequen¬ 
cies  in  the  range  of  3  to  60  Hz.  The  responses  of  critical 
structural  parameters  were  monitored  in  real  time  by  using 
a  computer-controlled  digital  structural  analyzer.  Addi¬ 
tionally,  the  effect  of  bandwidth  on  the  transfer  function 
was  studied  around  a  center  frequency  of  16  Hz.  It  was 
found  that  4  Hz  provided  an  adequate  bandwidth  for  the 
study  of  research  rotors.  The  current  S-61  rotor  has  a 
5/rev  frequency  of  16.9  Hz;  the  Sikorsky  UH-60  has  a 
4/rev  frequency  of  17.1  Hz,  and  the  Boeing-Vertol  360 
rotor  has  a  4/rev  frequency  of  17.6. 

The  shakers  were  controlled  by  a  digital  signal  gener¬ 
ator  whose  output  frequency  could  be  controlled  by 
microcomputer.  The  computer  was  programmed  to  con¬ 
duct  a  logarithmic  sweep  of  the  frequency  band  at  a  rate 
that  would  allow  the  RSRA  airframe  to  reach  steady-state 
response  at  any  given  frequency.  During  each  test,  the 
time-responses  of  all  transducers  were  recorded  on  mag¬ 
netic  tape,  using  the  Piloted  Aircraft  Data  System  (PADS) 
in  the  RSRA.  These  tapes  were  then  removed  from  the 
aircraft  and  processed  off-line,  as  described  below. 

The  time  records  of  each  transducer  and  the  master 
shaker  were  input  to  a  commercial  structural  analyzer, 
which  performed  a  fast  Fourier  transform  (FFT)  on 
the  signals  and  produced  a  least-squares  estimate  of 
the  complex  transfer  function.  This  transfer  func¬ 
tion  was  then  displayed  on  the  screen  as  a  plot  of 
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transduccr-rcsponsc/shaker-force  versus  frequency.  The 
data  to  create  the  plot,  in  the  form  of  256  real-imaginary 
pairs,  were  transferred  to  a  microcomputer  that  had  been 
programmed  to  interrogate  the  analyzer  and  record  the 
data  it  received  back  onto  a  digital  tape  cartridge.  Each 
transducer  produced  a  file  of  measurements  for  each 
degree  of  freedom  for  each  transducer  at  each  force  level 
tested. 

Using  commercially  available  software,  the  micro¬ 
computer  transmitted  the  reduced  data  to  a  mainframe 
computer,  for  processing  and  permanent  storage.  Main¬ 
frame  post-processing  was  divided  into  four  steps: 
(l)data  correction;  (2)  determination  of  slopes  and 
asymptotes  by  least-squares;  (3)  production  of  calibration 
matrices;  and  (4)  comparison  of  “predicted”  excitation  to 
actual  excitation. 

The  transfer  function  data  were  computed  relative  to 
twice  the  master  shaker  force.  The  slave  shaker  was  not, 
however,  perfeedy  in  phase  with  the  master  (a  lag  of  up  to 
about  6°),  nor  did  the  force  match  perfectly.  To  account 
for  this  effect,  the  real  and  imaginary  parts  of  these 
transfer  function  data  were  multiplied  by  the  complex 
factor  (Fm  +  Fs)/(2  x  Fm)  =  (1  +  Fs/Fm)/2.  To  extract 
slave-cylinder  data,  a  transfer  function  was  taken  of  the 
slave  cylinder  versus  the  master  cylinder  for  each  load 
level  and  excitation  axis  (shaker  inaccuracy  was  found  to 
be  a  function  of  both).  These  Fs/Fm  data  were  transmitted 
to  the  mainframe  computer  in  the  same  manner  as  the 
transducer  transfer-function  data.  A  mainframe  program 
then  read  the  Fs/Fm  data  and  applied  the  proper  correc¬ 
tions  to  the  corresponding  transducer  data. 

The  factors  co,  ci,  and  C2  from  equations  (13)  and 
(14)  were  determined  from  transformed  least-squares 
formulas  where 
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in  which  k  is  an  index  of  experiments  (different  force 
levels),  and  all  arithmetic  is  complex.  The  coefficients  co 
and  C2  for  each  transducer  (index  j)  and  degree  of  freedom 
(index  i)  were  computed  for  256  frequency  intervals  over 
the  range  of  14.5  Hz  to  18.5  Hz.  Each  coefficient  was  then 
weighted  by  normalizing  it  to  the  full-scale  value  of  the 
corresponding  transducer  to  prevent  ill-conditioning  of  the 
transfer  matrices. 

For  each  set  of  sensors  used  (i.e.,  for  each  case),  two 
calibration  matrices  were  formed,  one  using  co  as  the 
estimate  for  Hjj,  the  other  using  C2-  Each  such  matrix  was 
pseudoinverted  to  produce  a  candidate  calibration  matrix 
for  each  frequency  interval.  These  matrices  were  then 
multiplied  by  normalized  transducer  transfer  function  data 
to  produce  a  normalized  prediction  of  the  input  vector  (a 
vector  of  order  10  containing  real  and  imaginary  parts  for 
the  five  input  degrees  of  freedom).  This  was  then  com¬ 
pared  with  the  ideal  normalized  vectors  as  given  in 
table  2. 

An  “error  radius,”  £(&)),  was  defined  as  an  ad  hoc  measure 
of  goodness-of-fit  for  the  matrices.  It  is  defined  as  the 
square  root  of  the  sum  of  the  squares  of  all  errors  in  the 
10-dimensional  measurement  space: 
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TABLE  2.-  IDEAL  NORMALIZED  INPUT  VECTOR 


Input  DOF 

X 

Y 

z 

L 

~  M  | 

Re 

Im 

Re 

Im 

Re 

Im 

Re 

Im 

Re 

Im 

Longitudinal 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Lateral 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

Vertical 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

Roll 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

Pitch 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

Combined 

0.612 

0 

0.612 

0 

0.5 

0 

0 

0 

0 

0 
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where 


|F(0))| 


(20) 


The  term  error  radius  arises  because  its  value  repre¬ 
sents  the  magnitude  of  a  vector  from  the  actual  force 
vector  to  the  measured  force  vector  in  a  space  of  dimen¬ 
sion  10.  Since  the  single  DOF  data  were  used  to  make  the 
calibration  matrices  in  the  first  place,  the  error  radius  rep¬ 
resents  only  (1)  computational  errors,  such  as  roundoff 
error;  (2)  mathematical  uncertainty  resulting  from  ill- 
conditioning  and  stiffness;  and  (3)  nonlinear  behavior  of 
the  physical  system. 

Measurement  uncertainties  do  not  show  up  in  the 
error  radii  except  for  those  generated  from  the  multi-axis 
input  data,  which  were  not  used  to  make  the  calibration 
matrices.  Therefore,  differences  in  the  error  radii  between 
the  multi-axis  and  single-axis  data  sets  give  a  rough 
measure  of  the  effects  of  nonlinearities  in  the  structure 
versus  inaccuracies  in  the  transducer  data. 


RESULTS 

The  corrected  transfer  functions  of  the  individual 
transducers  at  different  excitation  levels  were  plotted  in  a 
three-dimensional  format  to  present  an  integrated  repre¬ 


sentation  of  the  response  of  each  transducer  to  shaker 
excitation.  Each  transfer  function  appears  on  the  plot  as  a 
wavy  or  jagged  line  connecting  the  256  data  points  in  the 
test-frequency  interval.  Each  transfer-function  curve  is 
connected  to  the  corresponding  point  in  the  neighboring 
curve  to  form  constant-frequency  curves  in  the  y-z  plane, 
giving  the  appearance  of  a  surface.  On  such  a  plot,  gross 
nonlinearities  show  up  very  prominently.  Force-dependent 
nonlinearities  (such  as  ankylosis)  would  appear  as  shown 
previously  in  figure  7,  and  frequency-dependent  nonlin¬ 
earities  (such  as  nonlinear  modes)  would  appear  as  shown 
in  figure  9. 

The  data  exhibited  the  expected  nonlinear  behavior 
described  above  in  many  cases,  as  may  be  seen  in  the  typ¬ 
ical  lest  data  shown  in  figure  13.  The  real  and  imaginary 
parts  of  measured  transfer-function  data  were  plotted 
against  the  input  shaker  force.  Curves  of  the  form 
(eq.  (14)) 


and(eq.  (15)) 


Hji  =  C2 

were  fitted  to  the  data  as  shown  above.  It  may  be  seen 
here  that  the  explanation  and  modeling  of  the  nonlinear 
behavior  given  above  is  at  least  reasonable  in  this  case. 
The  transfer  function  plots  for  the  28  transducers  are 
given  in  the  appendix.  They  are  grouped  by  force  direc¬ 
tion  and  ordered  by  transducer  as  in  table  1. 
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Figure  13  -  Asymptote  plot:  MRLIFTB  response  to  lateral  forcing. 


The  transfer  function  plots  for  the  load  cells 
(MRDRAG  through  MRGBQCF)  show  that  the  load  cells 
give  fairly  linear  responses  in  the  range  of  about 
17  through  18.5  Hz,  as  evidenced  by  the  fairly  smooth  and 
level  contours  of  the  plots  in  this  region.  At  the  low  end  of 
the  frequency  scale,  however,  a  force-dependent  nonlin¬ 
earity  is  seen  in  the  form  of  a  saddle-shaped  feature  cen¬ 
tered  on  about  200  lb  (890  N).  This  feature  is  seen  to  a 
greater  or  lesser  degree  in  all  the  load  cells  and  at  different 
values  of  the  transfer  function  (given  in  the  plots  as 
lb/1000  lb  of  excitation).  It  is,  therefore,  not  believed  to 
be  associated  with  the  load  cells  themselves,  but  with  the 
transmission  and  mounting  system. 

The  plots  for  the  transmission  accelerometers 
(LOMGB  through  XMRFPLO)  show  some  of  the  same 
features  as  those  for  the  load  cells.  The  saddle  structure 
noted  earlier  is  seen  in  these  plots  as  well,  also  centered 
on  200  lb,  though  in  some  cases  (XMRFBPLO  as  an 
example),  the  feature  is  inverted  to  form  a  hump.  As  gen¬ 
eral  observations,  note  the  following. 


1.  The  accelerometer  response  (shown  in  g’s/ 
1000  lb  of  excitation)  is  low,  probably  because  of  the 
large  mass  of  RSRA,  relative  to  the  excitation  levels 
allowable  for  this  test. 

2.  The  vertical  accelerometers  exhibit  very  good 
linearity  in  the  range  from  16.5  to  18.5  Hz. 

3.  Accelerometers  with  low  response  show  reduced 
linearity/repeatability.  Transducer  noise  becomes  signifi¬ 
cant  because  of  a  mismatch  (intentional,  to  simulate  flight 
hardware)  of  accelerometer  range  and  measured  accelera¬ 
tion.  Transducer  noise  shows  up  as  jaggedness  in  the  plot 
and  is  most  pronounced  for  low  force  inputs. 

4.  The  accelerometers  that  have  a  high  enough 
response  to  be  reliable  give  plots  with  striking  similarity 
to  load-cell  plots  for  load  cells  of  similar  orientation  (e.g., 
compare  XMRFBPV  and  MRLIFTC  or  compare  LOMGB 
and  MRDRAG). 
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The  results  support  the  hypothesis  that  the  main  gearbox 
and  frame  assembly  respond  as  a  rigid  body  relative  to  the 
response  of  the  fuselage  (as  explained  below). 

The  15  airframe  accelerometers  must  be  examined 
with  their  locations,  relative  to  the  point  of  excitation,  in 
mind.  Accelerometers  directly  below  the  transmission, 
XMRFRV  and  XMRFLV,  show  similar  responses  to  the 
transmission  load  cells  and  accelerometers.  Accelerome¬ 
ters  at  the  airframe  extremes,  such  as  VWGTPRT  or 
VWGTPLT  on  the  wings,  or  LATTPYLN  or 
VTTPYLON  atop  the  vertical  tail,  show  responses  domi¬ 
nated  by  airframe  structural  dynamics  and  bear  little 
resemblance  to  the  transmission  transducers.  Note  espe¬ 
cially,  the  structural  mode  at  17.5  Hz,  which  is  absent 
from  transmission  and  wing  responses.  Intermediate 
accelerometers,  such  as  S222FLRV  or  S450FLV,  show  a 
response  between  these  extremes.  These  observations 
support  the  expected  conclusion  that  the  fuselage  is  more 
flexible  and  behaves  less  like  a  rigid  body  than  does  the 
transmission. 

The  “local  modes”  seen  in  the  transfer  function  plots 
for  LATTPYLN  and  VTTPYLON  are  worthy  of  note. 
Modal  analysis  holds  that  all  modes  of  vibration  are 
global  properties  of  a  structure  and  exist  everywhere  in 
the  structure,  but  in  this  case  the  responses  of  all  the  other 
transducers  are  negligibly  small  and  the  mode  may  be 
thought  of  as  localized  in  the  tail  pylon.  Fabunmi  (ref.  7) 
has  shown  that  the  accuracy  of  the  force-determination 
method  is  limited  when  one  mode  dominates  the  response 
of  a  structure  (i.e.,  at  or  near  the  natural  frequency  of  the 
dominating  mode). 

Furthermore,  these  modes  exhibit  nonlinear  behavior. 
Note  that  the  constant-frequency  lines  near  the  classical 
mode  in  LATTPYLN  are  far  from  level;  they  rise  and  fall 
as  the  mode  shifts  in  frequency  with  increasing  force 
input  These  effects  are  difficult  to  model,  since  they  do 
not  affect  each  constant-frequency  cross  section  in  the 
same  manner,  and  are  not  accounted  for  at  all  using 
Dynamic  Force  Determination. 

Calculated  error  radii  for  cases  1  through  5,  for  use  of 
C2  coefficients,  are  shown  in  figures  14  through  18.  The 
individual  transducer  plots  are  presented  in  the  appendix. 
Case  5  (fig.  18)  represents  the  best  performance  of  the 
measurement  system  for  all  five  input  axes.  Each  data 
point  represents  the  magnitude  difference  between  pre¬ 
dicted  force  and  actual  force  normalized  to  100%  of  the 
nominal  applied  force.  Each  of  the  first  five  curves  for 
each  case  represents  the  accuracy  of  the  method  for  a  par¬ 
ticular  degree  of  freedom  and  reveals  only  how  well  the 
calibration  matrices  were  able  to  “predict”  the  statistical 


data  that  generated  them,  as  explained  above.  The  excep¬ 
tion  to  this  is  the  “combined,”  or  multi-axis,  load  data, 
which  represent  data  outside  this  statistical  universe. 


Figure  14.-  System  error  radius:  Case  1,  load  cells  only. 


Figure  15.-  System  error  radius:  Case  2,  load  cells  and 
gearbox  accelerometers. 
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Figure  16.-  System  error  radius:  Case  3,  load  cells  and 
airframe  accelerometers. 


Force,  lb 

Figure  17.-  System  error  radius:  Case  4,  airframe 
accelerometers  alone. 

The  predictions  are  least  accurate  at  low  force  levels, 
but  begin  to  level  out  at  200-300-lb  input  force  (lb-ft  for 
moments).  Here  the  ankylosis  effect  diminishes.  Beyond 
this  point  each  curve  has  a  minimum  in  the  400-  to  700-lb 
range.  This  is  thought  to  be  an  artifact  of  the  model  for 
obtaining  Hj;  coefficients,  and  not  to  represent  some 
“ideal”  force  level  for  the  system. 


Figure  18  -  System  error  radius:  Case  5,  all  transducers. 
CONCLUSIONS 

It  seems  clear  from  the  size  of  the  error  radii  that  a 
full-scale  dynamic  calibration  using  Dynamic  Force 
Determination  would  only  result  in  a  vibratory  measure¬ 
ment  capability  with  at  best  a  5%  accuracy.  This  limit  is 
imposed  by  the  RSRA  airframe  and  by  the  method  used. 
The  best  accuracy  achieved  with  data  not  used  to  form 
calibration  matrices  was  of  the  order  of  20%.  This  is  the 
best  accuracy  that  can  be  reasonably  expected  with  the 
configuration  tested  here.  This  is  far  less  than  the  accu¬ 
racy  expected  of  a  research  tool. 

A  dynamic  calibration  is  not  ruled  out  altogether,  if 
the  method  is  accurate  enough  to  reliably  reduce  the  error 
radii  to  under  5%  of  measured  force  or  if  a  different  mea¬ 
surement  or  analysis  technique  is  identified  and  used. 
There  is  hope  on  both  fronts. 

The  Dynamic  Force  Determination  used  here  was 
based  on  a  linear  model  of  airframe  response.  The  method 
may  be  extended  to  a  nonlinear  model  by  using  coeffi¬ 
cients  from  a  curve  fit  of  transducer  data  and  applying 
them  to  a  predictor-corrector  numerical  algorithm,  in 
which  the  linear  model  already  developed  becomes  the 
force  predictor  for  the  nonlinear  force  corrector.  Devel¬ 
opment  of  such  a  model  could  extend  the  accuracy  to 
lower  force  levels.  It  is  uncertain  at  this  time  whether  the 
approach  would  exhibit  enough  numerical  stability  to  be 
reliable  in  all  cases,  and  whether  the  increased  computer 
time  required  would  be  justified  by  the  results. 

Work  along  different  lines,  using  a  nonlinear  model 
of  load-cell  deflection  angles  and  parameter  estimation 
techniques  on  time-domain  data  has  been  done  and  was 
reported  by  others. 
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APPENDIX 


Graphs  that  represent  the  responses  of  28  transducers 
in  the  RSRA  to  excitation  in  five  degrees  of  freedom  at 
the  rotor  head  and  one  coupled  excitation  in  three  axes 
(multi -ax  is  mode)  are  presented  at  the  end  of  this  report. 
The  explanation  of  the  axes  in  these  “transducer- 
response”  plots  is  given  in  the  main  body  of  the  repon  in 


the  Results  section.  They  are  ordered  by  transducer  as 
shown  in  table  1  and  grouped  by  force  direction. 

The  author  gratefully  acknowledges  the  fine  work 
contributed  by  Cathie  M.  Jacobson  in  processing  the  data 
and  producing  the  graphs. 
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(b)  MRLIFTA  Concluded. 
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Figure  1.- Continued. 


Real  Imaginary 


Figure  1 Continued. 


ase 


QQQV/t)')  V»0\\0Vit\i 


Magnitude  Phase 


V*°c\ 


37 


(j)  VMGB  Concluded. 


se 


Magnitude  Phase 


(p)  VWGTPRT  Concluded. 


Rea!  Imaginary 


1 


^  UOVVSVKN*  i*ysuQiv 


OQi^V/^  viO\oviT\^  j,*\»uo*v 


50 


Figure  1.-  Continued. 
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(q)  VWGTPLT  Concluded. 


00<i 


<u 

</l 

o 

_c 

Q. 


41 

TJ 

3 

C 

cn 

O 

2 


<h  ‘»»ovva 


V.*<V\  ^*\*v»'o»v 


T3 

<D 

*o 

J3 

”o 

c 

o 

U 

z 

o 

J 

a, 


•8 

3 

C 

n 

O 


U 


e 

=3 

M 

£ 


59 


Magnitude  Phase 


(v)  XMRFLV  Concluded. 
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(w)  XMRFRV  Concluded. 
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Figure  1 Continued. 
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Figure  2.-  Continued. 
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Figure  3.-  Continued. 
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Figure  4.—  Continued. 
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(h)  LOMGB  Concluded. 
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